home *** CD-ROM | disk | FTP | other *** search
/ QRZ! Ham Radio 8 / QRZ Ham Radio Callsign Database - Volume 8.iso / pc / files / dsp / 96000tar.z / 96000tar / 96000 / appb / b119.asm < prev    next >
Assembly Source File  |  1992-04-28  |  4KB  |  89 lines
  1. ; This program was originally published in the Motorola DSP96002 Users Manual
  2. ; and is provided under a DISCLAIMER OF WARRANTY available from Motorola DSP
  3. ; Operation, 6501 William Cannon Drive West, Austin, Texas 78735-8598.  For
  4. ; more information, refer to the DSP96002 Users Manual, Appendix B, DSP
  5. ; Benchmarks.
  6. ;
  7. ; B.1.19    General Lattice Filter  
  8. ;                                                  GENERAL LATTICE 
  9. ; 
  10. ;              COEFFICIENT AND STATE VARIABLE STORAGE 
  11. ;       r0                         r4 
  12. ;                                    
  13. ;     
  14. ;    x: k3 k2 k1 w3 w2 w1 w0    y: s4 s3 s2 s1 
  15. ; 
  16. ;       m0=6 (=2*N, mod 7)     m4=3 (=N, mod 4) 
  17. ; 
  18. ;                                                SINGLE SECTION 
  19. ; 
  20. ;                                             EQUATIONS: 
  21. ;         
  22. ;                                          t'=t-k*s 
  23. ;                                          s'=s+k*t' 
  24. ;                                           t'?t 
  25.                       
  26. ;                                          output= sum(s'*w) 
  27. ;                  
  28. ;                        
  29. ;            
  30.           
  31.               
  32.                
  33.              
  34.  
  35. ;                                           DSP56000 IMPLEMENTATION 
  36. ;                                                                              Program    ICycles
  37. ;                                                            Words 
  38. ; move   #k,r0              ;point to coefficients 
  39. ; move   #2*N,m0            ;mod 2*(# of k's)+1 
  40. ; move   #state,r4          ;point to filter states 
  41. ; move   #N,m4              ;mod on filter states 
  42. ; movep  y:datin,a          ;get input sample 
  43. ; move          x:(r0)+,x0 y:(r4)-,y0 ;get first k, first s     1    1 
  44. ; do   #N,_el                         ;do filter                2    3 
  45. ; macr -x0,y0,a            b,y:(r4)+  ;t-k*s, save prev s       1    1 
  46. ; move          a,x1       y:(r4)+,b  ;copy t',get s again      1    1 
  47. ; macr  x1,x0,b x:(r0)+,x0 y:(r4)-,y0 ;t'*k+s,get k,get s       1    1 
  48. ;_el 
  49. ; move                     b,y:(r4)+  ;sv scnd to 1st st        1    1 
  50. ; clr   a                  a,y:(r4)+  ;save first state         1    1 
  51. ; move                     y:(r4)+,y0 ;get last state           1    1 
  52. ; rep   #N                                                      1    2 
  53. ; mac   x0,y0,a x:(r0)+,x0 y:(r4)+,y0 ;do fir taps              1    1 
  54. ; macr  x0,y0,a              (r4)+    ;finish, adj pointer      1    1 
  55. ; movep a,y:datout                    ;output sample 
  56. ;                                                             ---   --- 
  57. ;                                                                             Totals: 12  4N+10 
  58. ; 
  59. ;                                        DSP96002 IMPLEMENTATION 
  60. ;                                                            Program    ICycles
  61. ;                                                            Words 
  62. ; 
  63.  move   #k,r0              ;point to coefficients 
  64.  move   #2*N,m0            ;mod 2*(# of k's)+1 
  65.  move   #state,r4          ;point to filter states 
  66.  move   #N,m4              ;mod on filter states 
  67.  
  68.  move p  y:datin,d1         ;get input sample 
  69.  
  70.  move        #2,n4                                        ;  1       1 
  71.  move                        x:(r0)+,d5.s y:(r4)-,d6.s    ;  1       1 
  72.  do   #N,_elat                                            ;  2       3 
  73.  fmpy d5,d6,d0  fadd.s d0,d3                              ;  1       1 
  74.                 fadd.s d0,d1 d6.s,d3.s    d3.s,y:(r4)+n4  ;  1       1 
  75.  fmpy.s d5,d1,d0             x:(r0)+,d5.s y:(r4)-,d6.s    ;  1       1 
  76. _elat 
  77.                 fadd.s d0,d3                              ;  1       1 
  78.  fclr   d0                                d3.s,y:(r4)+    ;  1       1 
  79.  fclr   d1                                d1.s,y:(r4)+    ;  1       1 
  80.  move                                     y:(r4)+,d4.s    ;  1       1 
  81.  rep   #N                                                 ;  1       2 
  82.  fmpy d5,d4,d0  fadd.s d0,d1 x:(r0)+,d5.s y:(r4)+,d6.s    ;  1       1 
  83.                 fadd.s d2,d3              (r4)+           ;  1       1 
  84. ; 
  85.  move p  d3.s,y:datout         ;output sample 
  86. ;                                                           ---     --- 
  87. ;                                                   Totals:  14    4N+12 
  88.  
  89.